Spin Transfer Torque with Spin Diffusion in Magnetic Tunnel Junctions 
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Spin transport in magnetic tunnel junctions in the presence of spin diffusion is considered theo- 
retically. Combining ballistic tunneling across the barrier and diffusive transport in the electrodes, 
we solve the spin dynamics equation in the metallic layers. We show that spin diffusion mixes the 
transverse spin current components and dramatically modifies the bias dependence of the effective 
spin transfer torque. This leads to a significant linear bias dependence of the out-of-plane torque, 
as well as a non-conventional thickness dependence of both spin torque components. 

PACS numbers: 75.60. Jk,85.75.Dd,72.25.-b 
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Current-driven control of the magnetization direction 
in magnetic nanodevices has been made possible by 
the prediction of Spin Transfer Torque (STT), by Slon- 
czewski and Berger— . The transfer of spin angular mo- 
mentum between the spin-polarized electrical current and 
the local magnetization has been observed in various de- 
vices such as metallic spin-valves, magnetic tunnel junc- 
tions (MTJs) and magnetic domain walls* 2 * 3 -. A number 
of technological applications have been proposed, based 
on these devices, such as STT-MRAM, Race track mem- 
ory etcA 

The most promising candidate for memory applica- 
tions to date is the MTJ 5 due to its high perfor- 
mances and good compatibility with the existing C-MOS 
technology 6 . Therefore, understanding the nature of the 
spin transfer torque in such devices is of seminal im- 
portance. The first theoretical studies showed that in 
the ballistic limit of symmetric MTJs comprising semi- 
infinite ferromagnets, the spin torque is on the former— 

T = T||M x (P x M) + T ± M x P, (1) 
T,| = ai V + a 2 V 2 , T x = b + b 2 V 2 . (2) 

The first term Ty in Eq. (Q}, referred to as the in-plane 
torque, competes with the magnetic damping allowing 
for self-sustained magnetic precessions and switching, 
whereas the second term T±, referred to as the out-of- 
plane torque, acts like an effective field applied along 
P. Here, M and P are the magnetization direction of 
the free and pinned layers, respectively and V is the 
bias voltage applied across the junction. Interestingly, 
as shown in Eq. ([2]), the bias dependence of the spin 
torque components is well defined, the in-plane torque 
possessing both linear (eti) and (small) quadratic com- 
ponents (02 *C &r), and the out-of-plane torque being 
only quadratic (62), besides the zero-bias exchange cou- 
pling term (bo). This bias dependence has been con- 
firmed experimentally in FeCoB/MgO/FeCoB by spin- 
diode measurements^. The quadratic bias dependence 
of the out-of-plane torque has serious implications on 



current-driven magnetization dynamics. Whereas the 
sign of the in-plane torque depends on the polarity of 
current injection, the out-of-plane torque is always in the 
same direction. Consequently, at positive polarity both 
torques favor the antiparallel configuration, whereas at 
negative polarity, the in-plane torque favors the parallel 
configuration while the out-of-plane torque favors the an- 
tiparallel one. This competition leads to back-hopping of 
the magnetization state which is detrimental for applica- 
tions such as MRAM±^. 

However, recent experiments have reported important 
discrepancies between the actual bias dependence of the 
spin torque and the one proposed in Eqs. fl])-© 1 - - — 
In particular, Oh et al^ showed that in an asymmetri- 
cally designed MTJ, the bias dependence of the out-of- 
plane torque acquires a linear contribution. This con- 
tribution reduces the competition between in-plane and 
out-of-plane torques for negative polarity and results in a 
net reduction of the back-hopping process. Several mech- 
anisms have been proposed to alter the bias dependence 
of the spin torque, such as finite layer thicknes s 10 ' 11 ' 13 , 
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barrier and electrode asymmetry 
and spin-flip scattering in the electrode s 19 ' 20 . Neverthe- 
less, these mechanisms fail to explain the large linear bias 
dependence usually observed in the out-of-plane torque. 

In the present study, we propose a model of spin trans- 
fer in a diffusive spin transport approach taking place 
in the metallic electrodes of MTJs. We will show that 
the diffusion processes (i) mix the two spin components 
transverse to the local magnetization which has the result 
to (ii) dramatically modify the bias dependence of the ef- 
fective spin torque as well as its thickness dependence. In 
particular, we show that the out-of-plane torque acquires 
a significant linear bias dependent term. This result is 
obtained by solving the spin diffusion equation for the 
spin accumulation transverse vector in the metallic lay- 
ers adjacent to the insulator and imposing the interfa- 
cial spin current as a boundary condition. In particu- 
lar, our model adopts the point of view of Zhang et al. 
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whereby magnon- assisted tunneling process, at high bias, 
quenches any transport of hot electrons over a short dis- 
tance inside the ferromagnetic electrode resulting in a 
pure diffusive transport within the electrode^. 

Let us consider a MTJ composed of N/F/I/F/N stack, 
where I is the insulating spacer, F are the ferromagnetic 
layers and N the non magnetic metallic leads connected 
to the F layer. The tunneling process through the insulat- 
ing barrier imposes a ballistic injection of carriers at I/F 
interface which can be described by several calculation 
techniques like the free electron model, a tight-binding 
treatment or a density functional theory. These different 
approaches are well documented in a recent review paper 
Ref. [23l From the above arguments, it results that the 
two components of the transverse spin current Jo are im- 
posed at the I/F interface. It constitutes a boundary con- 
dition to the coupled diffusive spin transport-relaxation 
equations for the transverse component of the spin ac- 
cumulation vector m. Along these guidelines, the spin 
dynamics of the transverse spin accumulation in the fer- 
romagnetic electrode from the I/F interface is governed 
by the following time-dependent coupled equations 

Q -1-1 

= -V • J m x M M x (m x M) - — h) 

Ot Tj T T sf 



J = -VV <g> m, 



(4) 



where m is the spin accumulation, M is the direction of 
the localized 3d magnetization, t is the time, and J is 
the spin current tensor. In Eq. @, the expression of 
the spin current is limited to Ohm's law via the diffusion 
constant T>. These equations account for spatial variation 
of spin current, spin precession, spin dephasing and spin 
relaxation through the respective spin precession time 
tj, spin decoherence time and spin relaxation time 
r s f. Whereas the spin relaxation affects the three spin 
components, the spin precession and spin decoherence 
terms only affect the two transverse components of the 
spin accumulation vector. In the transient regime of spin 
injection dynamics, these equations can be solved analyt- 
ically using Green's functions techniques in presence of 
a source term JoS(z)5(t) (5 is the Dirac distribution and 
z is the direction of the tunneling current where z = 
correspond to the exact I/F interface position). In the 
steady-state regime of spin injection (x J Q dt), it simply 
gives 

V • J = — -m x M - — M x (m x M) — m, (5) 

tj 70 t s f 

J = -DV®m, (6) 

The solutions for the spin accumulation in the ferromag- 
net, when the magnetization is along z reads 

m t = my + im± = Ae x/L + Be~ x/L 
1 _ i 1 1 
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(7) 
(8) 



where Af = XVj. We are now interested in the spin ac- 
cumulation in the ferromagnetic layer, displayed in Fig. 



[TJ On the left interface I/F, the boundary conditions are 
given by the interfacial spin current injected through the 
barrier Jq = J® + ij±- The two components are the 
in-plane and out-of-plane tunneling spin currents defined 
by Jo = jftM x (P x M) + J°M x P. On the right 
interface F/N, the connection is given by the continuity 
of the spin accumulation and spin current. In principle, 
interfacial resistance and spin-flip give additional contri- 
butions to the interfacial boundary conditions. In the 
present model, we disregard these effects since they can 
be added by simply inserting an effective interfacial layer. 
After some algebra, the two components of the transverse 
spin density read 

j oLs inh^+7?cosh^ x ^ d 
T>f cosh y + 77 sinh ^ 
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JoL rje w 

- vm± = — -, -r, x > a 

L> F cosh -| + 77 sinh ^ 



(10) 



where rj = -jt- , T>i is the conductivity of the i-th layer 
and d is the thickness of the ferromagnetic layer. Figure 
[ljdisplays the spatial profile of the spin accumulation rn\\ 
and m_L transverse to the local magnetization for differ- 
ent spin dephasing length A^=0.5 nm, 1 nm, and 1.5 nm. 
In this figure, we assumed J^_=Q (no ballistic out-of-plane 
spin current). For the chosen parameters, the spin accu- 
mulation is not fully absorbed in the ferromagnet and 
decays away from the I/F interface. Furthermore, due to 
the spin precession, the spin accumulation displays both 
in-plane and out-of-plane components. 
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FIG. 1. (Color online) Transverse spin accumulation as 
a function of the distance in the electrode in F(2 nm)/N. 
The parameters are Aj=l nm,Af^ = 15 nm, nm and 

V F /V N = 1. 

Note that the out-of-plane component m± increases 
when enhancing the spin dephasing length. Therefore, 
even in the absence of out-of-plane interfacial spin cur- 
rent (j7^=0), the interplay between spin dynamics and 
spin diffusion is expected to produce a torque with both 
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in-plane and out-of-plane components. The spin torque 
is denned as the spatial change of spin current, compen- 
sated by the spin relaxation term 



T = - / dn f-V-J-—m 

n Jn V T *f 



(11) 



Here, f2 is the volume of the magnetic layer. Therefore, 
the total spin torque exerted on the ferromagnet is 



Jo L 2 cosh j 



77 sinh ■£ — 1 



d L 2 



cosh ^ 



77 sinh 



(12) 



where -L = -3V + i. 

In the limit of infinite ferromagnetic layer thickness, 
the torque in Eq. (fT2"]) reduces to 



T|| = ^(((l + x)/? 2 +X 2 )^ + X 2 /3J? 



(13) 



T^ = ^(((l + x)P 2 + X 2 )Jl-X 2 PJ$), (14) 

where x = ^/r s f, /3 = tj/t,/ and £ = (1 + x) 2 /3 2 + x 2 - 
In the case of very short spin dephasing % <C 1, as in 
Fc/MgO/Fei 2 -, the spin torque reduces to the ballistic 
limit, T|| = JjP/d and Tl = J^/d. In the limit of infinite 
spin dephasing x ^> 1, such as the one considered in 
Ref. [l9| and |20L the spin torque arises from a mixture 
of J° and Jl, T,| = (J* + /3j£)/d(l + /3 2 ) and T ± = 
G7° - /3^P)/d(l + /3 2 ). Therefore, in the case of semi- 
infinite ferromagnets, a linear bias dependence of the out- 
of-plane torque arises only if both the spin dephasing and 
the spin precession are comparable to the spin diffusion 
length. 

The case of a finite ferromagnetic layer presents signifi- 
cant differences due to the adjacent normal metal. When 
the spin precession and dephasing occur on a length much 
smaller than the thickness of the ferromagnet, the trans- 
verse spin current is mostly absorbed at the interface be- 
tween the insulator and the ferromagnet and the spin 
torque bias dependence is described by Eqs. ((U)-©. 
Conversely, when the spin dephasing and precession ex- 
tend on a length comparable to the layer thickness, the 
transverse spin current is no more confined to the inter- 
face and extends on the layer thickness to the second in- 
terface with the normal metal, as shown in Fig. [TJ This 
new dynamics mixes the two transverse components of 
the spin current and is therefore responsible for the de- 
viation from the ballistic bias dependence shown in Eq. 
@ . Figure [5] displays the two components of the spin 
torque, T11 and T±, obtained from Eq. (|12p in the case 
J7°=0 as a function of the thickness of the ferromag- 
netic layer. In this figure, the two contributions have 
been multiplied by the distance d to remove the 1/d de- 
pendence coming from the interfacial nature of the spin 
torque. Interestingly, it appears that both in-plane and 
out-of-plane torques possess a thickness-dependence that 
is not simply oc 1 /d, as it is expected in the case of purely 
interfacial spin torque. 
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FIG. 2. (Color online) Thickness dependence of the normal- 
ized spin torques dT« and dT± for different spin precession 
lengths, Xj — 0.5 nm, 1 nm and 1.5 nm. The parameters are 
the same as in Fig. [1] except for the spin dephasing length: 
Aj, — 0.5 nm (solid lines) and X^ — l nm (dotted lines). In- 
set: Thickness dependence of the ratio between out-of-plane 
torque and in-plane torque for the same parameters. 



A deviation from the 1/d-thickness dependence of the 
spin torque has been identified previously in calculations 
assuming a ballistic regim o 11 ' 12 . For example, Ref. [IH 
observes a deviation that is attributed to resonant quan- 
tum states in the finite ferromagnetic layer. Note that 
the model proposed in the present work assumes a diffu- 
sive transport in the metallic layers, i.e. defects and dis- 
order are strong enough to quench quantum coherence. 
Therefore, the deviation we observe here is not related 
to quantum states, but rather to the incomplete absorp- 
tion of the spin current: when the thickness of the free 
layer is on the order of or smaller than the spin dephas- 
ing length A^, the transverse spin current responsible for 
the spin torque is not fully absorbed in the free layer and 
diffuses towards the capping layer. This induces a devia- 
tion from the usual 1 / d-thickness dependence. Increasing 
the thickness of the free layer improves the absorption of 
the spin current and for thicknesses much larger than the 
spin dephasing length, the thickness dependence of the 
torque recovers the 1/d limit (see Fig. [5]). Note that 
in the case of half-metallic behavior, as in Fe/MgO/Fe 
tunnel junctions, the minority band with Ai symmetry 
does not propagate in the ferromagnet which results in 
a quenching of the spin dephasing length (d ^> A^) and 
reduces the spin torque to a 1/d behavior. 

Since the transverse spin accumulation is not fully ab- 
sorbed in the ferromagnet, the remaining unabsorbed 
transverse spin accumulation diffuses into the normal 
metal and modifies the actual spin torque exerted on 
the ferromagnet. More specifically, the mixing that gives 
rise to both components increases for smaller thicknesses. 
Therefore, the ratio between the out-of-plane and the in- 
plane components, Tj_/T||, increases when decreasing the 
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thickness of the ferromagnet, as displayed in the inset 
in Fig. [5] The nature of the spin torque asymptoti- 
cally tends towards the semi-infinite ferromagnet limit for 
thicknesses much larger than the spin dephasing length. 

Let us now illustrate the influence of the normal metal 
on the spin diffusion in the ferromagnetic layer. The role 
of the normal metal can be qualitatively understood if 
one models the F/N bilayer by an equivalent semi- infinite 
F layer with an effective spin diffusion length A*^ . In the 
semi-infinite limit, the spin torque reduces to 



*2 



Tn 



where 



.„ d->oo Jo L 



a:, 



(15) 



, x x ._. . „ s * being the effective 

J <p sf 

spin diffusion length arising from the presence of the ad- 
jacent normal metal. By equating Eq. (TT5T) with Eq. 
(d2), one can define the effective spin diffusion length 
A*^ as a function of the parameters of the normal layer. 
To obtain a tractable analytical result, we consider an 
ultrathin ferromagnet, so that 



.„ d-yO^ Jo L 
ll± > 



r2 , L 



This provides the analytical expression 
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where pjq is the probability to find the particle in N. 
The bulk spin relaxation time in the finite ferromagnet 
rfj is renormalized by the presence of the normal metal 
through the probability p^. Increasing the spin diffu- 
sion length of the normal metal A^- or decreasing the 
thickness of the ferromagnet d reduces the spin diffusion 
length in the ferromagnet. The second term of Eq. (TIE)) 
can be viewed as the effective spin lifetime in the finite 
ferromagnetic layer. 

This thickness dependence has a very important im- 
plication on the bias dependence of the spin torque. If 
one assumes a bias dependence of the spin current on 
the form J® = a\V and J® = b 2 V 2 , as expected and 

observed for systems such as Fe/MgO/Fe 12,14 , then both 
in-plane and out-of plane spin torque components will be 
a mixture of linear and quadratic bias dependences. For 
example, Sankey et al. and Kubota et al. experimentally 
found that T±/T\\ \200mV ~ 10% to 15% 14 , Assuming that 



the bias dependence of the interfacial spin current Jq is 
given by the spin torque bias dependence measured by 
Sankey et al. , the in-plane and out-of-plane torques are 
reproduced in Fig. [3] for different thicknesses of the ferro- 
magnetic layer. As expected from the discussion above, 
the linear character increases when decreasing the fer- 
romagnetic thickness and at large thicknesses it reduces 
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FIG. 3. (Color online) Bias dependence of the in-plane (a) 
and out-of-plane torques (b) for different thicknesses of the 
ferromagnetic layer. The parameters are Aj = A^ =1 nm, 
\ff=W nm, \^f=l nm and T>f /Cjv = 1- The bias dependence 
of the interfacial spin current j7o is taken from the spin torque 
bias dependence measured by Sankey et a/.— . The inset dis- 
play the zoom in the range ±200 mV. 



to the bulk value given by Eq. (|14l) . Since the ballistic 
out-of-plane spin current is small compared to the 
in-plane spin current JT^, only slight modification of the 
in-plane torque is expected. 

In conclusion, the influence of spin diffusion in the 
metallic layers of MTJs on the spin transfer torque has 
been addressed theoretically. Assuming an interfacial 
bias-driven spin current at the interface between the in- 
sulator and the ferromagnet, the spin diffusion equation 
is solved and describes a complex spin dynamics in the 
metallic layers. It is found that this dynamics mixes 
the components of the spin current tranverse to the lo- 
cal magnetization which results in a superposition be- 
tween linear and quadratic bias dependence for both the 
in-plane and out-of-plane torques. The thickness depen- 
dence of the spin transfer torque is also altered for small 
thicknesses. 

The authors acknowledge fruitful discussions with A. 
Fert. 
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